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From neutron inelastic and elastic scattering we have determined the magnetic structure and 
fluctuations in the YBa2Cu3 06.35 superconductor (Tc=18 K), near the boundary of the supercon- 
ducting phase. The long-range ordered coUinear spins of the insulating antiferromagnet are replaced 
by a commensurate central mode arising from slow, isotropically polarized, short-range spin corre- 
lations extending over four planar unit cells. The inelastic spectrum up to 30 meV is also broad in 
wave vector and commensurate. In contrast to the the resonance peak of higher T^ superconductors, 
the spins exhibit a single overdamped spectrum whose rate of relaxation F decreases on cooling and 
saturates at 2F=5 ± 1 meV below ~ 50 K. As the relaxation rate saturates the quasi-static spin 
correlations grow and become resolution limited in energy. The spin susceptibility above ~ 50 K 
follows the same lu/T scaling relation found for the monolayer La2_2,Sr3;Cu04 system indicating that 
the dominant energy scale is set by the temperature. Below 50 K the scale length is geometric and 
not linked by velocity to dynamic widths. Despite the large differences from an antiferromagnet, we 
show that integrated intensity conserves the total moment sum rule, and that on cooling the spectral 
weight transfers from the inelastic spin relaxation to the elastic central peak. There is no observ- 
able suppression of the spin fluctuations or central mode upon the onset of superconductivity. The 
spins respond not to coherent charge pairs but to hole doping allowing coexistence of glassy short 
range spin order with superconductivity. Since the physics of the weakly superconducting system 
YBC06.35 must connect continuously with that in more strongly superconducting YBCOe.s, we find 
that neither incommensurate stripe-like spin modulations nor a well-defined neutron spin resonance 
are essential for the onset with doping of pairing in a high temperature cuprate superconductor. 

PACS numbers: 74.72.-h, 75.25.+Z, 75.40.Gb 



I. INTRODUCTION 

There is a direct connection between high temper- 
ature superconductivity and antiferromagnetisniFii^i^i^ 
Despite much effort, however, a universal theory behind 
the cuprate phase diagram has not been estabhshed. One 
of the most interesting regions of the phase diagram lies 
near the antiferromagnetic phase and the onset of su- 
perconductivity. This region with its low carrier con- 
centration and highly suppressed transition temperature 
should provide much needed information on how antifer- 
romagnetism is destroyed in favor of high-temperature 
superconductivity. 

When the superconducting order parameter is larger, 
it is known from an extensive study (Ref. 00) in oxygen 
and structurally ordered YBa2Cu306.5 (YBCOe.s with 
Tc=59 K, corresponding to a hole doping of p=0.09), 
that the spin spectrum below and above Tc extends to 
the lowest energies but with no elastic magnetic Bragg 
peak. The results do not support expectations that a 
new hidden order parameter might exist that would have 
created a spin gap and associated Bragg peak. The low- 
energy spin response is modulated incommensurately in 
momentum and rises to form an hour-glass shaped spec- 
trum around a commensurate resonance at energy E=33 
meV. For temperatures below Tc the low-energy suscep- 



tibility is suppressed and its spectral weight is trans- 
ferred to the resonance energy thereby conserving the 
total moment sum rule. Above Tc the substantial broad 
response near resonance energies suggests incoherent pair 
formationi^ Experiments show that the resonance de- 
fines a crossover point that separates the low-energy in- 
commensurate region from the high-energy region where 
the dynamics resemble spin-wave excitations.^'^'^ Sim- 
ilar high-energy structure has been observed at larger 
doping in YBCOg-i-a; (Hayden et alr^) and in monolayer 
La2-2;Baj;Cu04 (Tranquada et al}^) and La2-a:Sra:Cu04 
(Christensen et al}^). Other studies of the spin excita- 
tion spectrum for similar hole doping in YUCOq+x by 
Fong et al}^, Dai et aL^'' and Pailhes et aL~ have mostly 
focussed on energies around the resonance. 

Studies of La2-a:Sra;Cu04 (LSCO) for low carrier dop- 
ing show that a spin-glass phase lies interposed between 
the antiferromagnetic insulating and the superconduct- 
ing phasesii^ii^ Such an intermediate phase has not been 
identified in YBCO. In LSCO the lightly doped antifer- 
romagnetic insulator was found to break up into islands 
of correlated spins with moments pointing within the 
Cu02 planes and the intensity from the ordered Cu^"*" 
moments decays rapidly on doping so as to enter the 
superconducting phase j^^'^^ This same phase and mag- 
netic structure has recently been confirmed for overdoped 
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La2-2;Sr2;Cu04+j, The superconducting and spin glass 
phases of LSCO both exhibit incommensurate peaks that 
are dynamic and static respectively. Spin coupling to 
the lattice is more important for LSCO than for YBCO 
in view of the abrupt strongly first-order transition with 
doping from diagonal to coUincar incommensurate wave 
vectors The spin response follows lu/T scaling relation 
with temperature, T, over a broad range of energy oj with 
a breakdown observed at low-temperatures from the out- 
of-plane anisotropyjSS, Inelastic scattering on Li doped 
La2Cu04^^'^^, where the Li causes a large perturbation 
as it replaces the Cu spin and charge in the supercon- 
ducting Cu02 plane, also shows a range of uj/T scaling. 
Surprisingly it was suggested that the inelastic proper- 
ties indicated the close proximity of a quantum critical 
point, in contrast to the spin-glass phase in LSCO and 
proximity to a known first order transition in the mag- 
netic structure. 

For YBCO there may be two classes of magnetic be- 
havior as Tc is suppressed. All studies for the fairly well 
doped YBCOg+x cuprates described above have shown 
that the magnetic spectrum differs considerably from 
that of the insulator. Its main features are the absence of 
elastic scattering, dynamic incommensurate fluctuations, 
a commensurate resonance, cones of highly damped spin 
waves and metallic resistivityi^'^'^ i^^'^'^i^^'^"^ It is therefore 
of great importance to reduce further the hole doping to 
where the resistivity is insulating and then to observe 
how the antiferromagnetic phase evolves from the super- 
conducting state and how the fluctuations are altered 
by the presence of mobile holes. In particular, does the 
low-energy spin spectrum begin to resemble that of the 
antiferromagnet as the doping falls to the critical value 
Pc=0.055 where superconductivity disappears? We will 
show that it does not in general. The similarities to the 
antiferromagnet are that the scattering becomes com- 
mensurate, that an elastic signal appears on cooling, and 
that the inelastic signal drops from a resonance region 
to very low energies. The differences are that the elastic 
and inelastic signals are broad in momentum signifying 
short-range correlations, that the elastic signal develops 
continuously without a transition, and that the spin di- 
rections remain isotropic. 

Remarkably little neutron inelastic scattering has been 
done in YBa2Cu306+2: (YBCOg+i) for charge dop- 
ing near the onset of superconductivity. In YBCOg.as 
(Tc=39 K) Mook et al^ observed incommensurate scat- 
tering at low-energies and a resonance at ~ 23 meV. Al- 
though the sample studied by Mook et al. had the same 
nominal oxygen concentration as the crystal we have in- 
vestigated, the in-plane doping was clearly much larger 
as shown by its large Tc=39 K. Its properties place it 
in the fairly well-doped class described above. The sam- 
ple we have studied has a more highly suppressed super- 
conducting transition temperature of Tc=18 K or ~0.2 
Tc{optimal) — 93K and a lower toin-plane hole doping 
and thus lies much closer to the boundary of the super- 
conducting phase. It belongs in the low-doped class. 



Most phase diagrams show that a cooled sample 
traverses a single boundary into the superconducting 
phase. However, data from spin resonance of the pos- 
itive muon26'27 (^SR) in YBCOe+x, in Ca doped YBCO 
and in LSCO, have been used to suggest a very different 
phase diagram. Two additional transitions are inferred 
within the superconducting dome, one to Neel order and 
at lower temperature a transition to a frozen glass. We 
will show using neutron elastic and inelastic scattering 
that in YBCOg.ss there is no transition to Neel order 
and that the nearly static Cu^"*" spins are isotropically 
polarized and have short correlation lengths. 

Transport measurements have recently revealed some 
very interesting properties around the critical hole con- 
centration for superconductivity. Sutherland et al. 
looked in YBCO and in LSCO for the electronic (linear- 
in-T) part of the thermal conductivity When they 
reduced the doping to leave the superconducting phase 
they found YBCO to remain metallic whereas LSCO be- 
came insulating. This marks a clear difference between 
the two systems over and above the spin coupling to ro- 
tation of the LSCO octahedra. Liang et al. observed 
a nonlinear relation between the lower critical field Hd, 
Tc and the superfluid density in extremely underdoped 
YBCOe+x, in contrast to expectations from BCS the- 
oryf22 Photoemission results also show a direct correla- 
tion between electronic properties and hole doping, some 
of which have been correlated with the spin fluctua- 
tions. Magnetoresistance measurements have found a 
large anisotropy in heavily underdoped YBCOg+i: that 
has been attributed, without observational evidence, to 
the formation of charge stripes. '^^ It is therefore impor- 
tant to carry out a detailed neutron study in this region 
of low hole doping to determine if there is a Cu^+ spin 
origin to the anomalous electronic properties. 

A general picture of how the spin structure and dynam- 
ics evolves as carrier concentration is reduced is given in 
Fig. [T]which illustrates the qualitative change in the neu- 
tron scattering between Ortho-II YBCOe.s (Tc=59 K) 
and YBCOg.as (Tc=18 K) well below their superconduct- 
ing temperatures. The spectral weight of the intensity, 
proportional to S(Q, uj) undergoes a dramatic change as 
doping is reduced. In YBCOe.s the response contains no 
elastic peak, and consists solely of a continuum termi- 
nating at a well-defined resonance energy of ~ 33 meV, 
while the response in YBCOe.ss has been replaced by 
an intense elastic peak and, at very low energy, a weak 
inelastic feature at ~ 2.5 nieV from spins that are re- 
laxing not resonating. We will show a direct connection 
between the two, in which the relaxation rate of the in- 
elastic feature slows on cooling and drives up the intensity 
of the central elastic component. There is no detectable 
incommensurability. As we shall see later, the total spec- 
tral weight, integrated over momentum and energy, is the 
same within error for the same energy range as that in 
the Ortho-II YBCOg.s superconductor, and therefore is 
the same as that of the antiferromagnet as shown earlier 
Thus the spin spectral weight is conserved as carriers are 



reduced towards the critical doping, despite the dramatic 
evolution of the spectral form. 

This paper provides an overview in energy and temper- 
ature of the spin response measured with thermal and 
polarized neutrons, and it provides the foundation on 
which part of a previous short paper was based. '^^ There 
we focused on the scattering only around huj—O. Here we 
fully describe the relationship between the inelastic spin 
relaxation and the central peak as a function of temper- 
ature. The first section gives the method, the second 
a detailed analysis of the central component, and the 
third section describes the higher energy spectrum up to 
energy transfers of 35 meV. In the final section we dis- 
cuss the integrated intensity and shows a direct connec- 
tion between the scattering centered around huj=0 and 
the broad inelastic scattering extending up to 30 meV at 
room temperature. 



II. EXPERIMENT 

Polarized and unpolarized neutron scattering mea- 
surements were made with the C5 DUALSPEC spec- 
trometer at the NRU reactor. Chalk River Laborato- 
ries. For unpolarized measurements a vertically focussing 
graphite (002) monochromator and a graphite (002) an- 
alyzer selected incident and scattered energies. A py- 
rolytic graphite filter in the scattered beam eliminated 
higher order reflections, and a liquid nitrogen cooled sap- 
phire filter before the monochromator reduced the fast 
neutron background. Elastic measurements were made 
with 2 graphite filters after the sample while inelastic 
measurements were made with 1 filter. Where we have 
compared elastic and inelastic intensities, we have cor- 
rected for the extra measured attenuation of the second 
filter. For unpolarized measurements with fixed E/=14.6 
meV the horizontal coUimation was [33' 29' S 51' 120']. 
For unpolarized measurements with E/=5 mcV the col- 
limation was [33' 48' S 51' 120'] and a liquid nitrogen 
cooled beryllium filter placed before the monochromator 
eliminated higher order neutrons. We have corrected all 
intensity the data for contamination of the incident beam 
monitor by higher wavelength neutrons as described else- 
where. ^''^^ 

To determine the magnetic structure and to extract the 
magnetic scattering over a broad range in energy trans- 
fer, we relied on polarized neutrons.^- The neutron beams 
were polarized with Heusler (111) crystals as monochro- 
mator and analyzer. A graphite filter suppressed higher- 
order scattered neutrons. A Mezei flipper in the incident 
beam allowed spin-flip (SF) and non-spin-flip (NSF) cross 
sections to be measured. Two pairs of coils applied to the 
sample a weak 3-5 Gauss field to control the neutron spin 
direction, either a horizontal field (HF) parallel to Q, or a 
vertical field (VF) perpendicular to Q. The flipping ratio 
was measured to be 12:1 for both configurations. With 
the flipper on, the difference between the HF and VF 
count rates gave the spin-flip scattering from the mag- 
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FIG. 1: The inelastic spectrum S(Q,tj) of YBCOe.as (Tc=18 
K) compared with that of Ortho-II YBCOe.s (Tc=59 K) 
both measured in the superconducting state. The momentum 
transfer is Q=(0.5,0.5). The spectral weight in YBCOe.ss is 
concentrated at much lower energies relative to YBCOe.s, but 
its integral is approximately conserved. Where no error in in- 
tensity is visible it is smaller than the symbol size because it 
is determined accurately from scans at constant energy. The 
closed (open) symbols are polarized (unpolarized) absolute 
intensities and the diamonds show the central mode intensity 
measured with polarized neutrons reduced by 3. 



netic electrons alone independent of nuclear incoherent 
or phonon scattering. To change field and prevent any 
depolarization from trapped flux, the sample was always 
warmed to ^ 50 K (more than twice Tc) and then field 
cooled to low temperatures in the new field. Both field 
coils (and hence the field direction) were fixed with re- 
spect to the sample. 

The sample consisted of seven ^ 1 cc crystals, mutu- 
ally aligned in the (HHL) scattering plane with an over- 
all mosaic of ~ 1.5°. From DC magnetization the critical 
temperature was measured to be Tc = 18 K with a width 
of ~ 2 K corresponding to a hole doping of p=0.06 (based 
on the formula provided by Tallon et aL— ). This dop- 
ing is close to that from a detailed study of hole doping 
and lattice constants by Liang et al.^^ The sample is or- 
thorhombic with lattice constants a=3.843, 6=3.871, and 
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FIG. 2: In plane elastic scattering scans through the (1/2, 
1/2, 2) magnetic peak at 2 K and 25 K with E/ = 14.6 meV. 
The magnetic intensity was obtained by subtracting a flat 80 
K background scan. The solid lines are results of fits to a 
Lorentzian squared as described in the text. The horizon- 
tal bar represents the resolution width. The inverse intrinsic 
width (and hence correlation length) does not diverge at low 
temperatures. 



c=11.788 A and is twinned. The same crystals, grown by 
a technique described previouslj*^, were used to obtain 
the low-energy properties with cold neutron and the high- 
energy spectrum with spallation neutrons. "^^'"'^ Neutron 
diffraction revealed no well-defined oxygen superlattice 
Bragg peaks indicating that any oxygen ordering in the 
chains is short ranged. 

To compare our data with those obtained from other 
oxygen concentrations and to determine the spectral 
weight, we normalized the intensity to an absolute 
scale by means of a constant-Q scan through a trans- 
verse acoustic phonon near the (006) Bragg peak at 
(0.15,0.15,6) at 85 K. This is the same procedure used 
previously where, by comparison to other calibration 
methods using a vanadium standard, we have found the 
systematic error in conversion to absolute units to be 
~ 15-20 %^ To confirm the cahbration via a phonon 
we have carried out a further independent calibration 
against the Bragg peak integrated intensity. This proce- 
dure, used elsewher o^^'^^ , is described in the Appendix. 
Having shown that the two methods agree well, we then 
quote the absolute intensities based on the phonon cali- 
bration. 



III. ELASTIC SCATTERING: CENTRAL PEAK 

We previously discovered that the inelastic spectrum 
at low temperatures exhibits two distinct energy scales, 



a slow central peak, energy-resolution limited as T— >0, 
and a fast timescale corresponding to broad overdamped 
inelastic fluctuations. We will show that the cen- 
tral peak exhibits three-dimensional (3D) correlations 
whereas the fast inelastic response displays only dif- 
fuse scattering from independent bilayers characteristic 
of two-dimensional (2D) correlations. The onset of the 
central peak can then be seen as a crossover from 2D 
dynamic correlations to quasi-static 3D correlations. 



A. Unpolarized Experiments 

With unpolarized neutrons we have determined the 
elastic magnetic scattering and in-plane correlations in 
scans of the form (H,H,2) and (H,H,5), and also deter- 
mined the interplanar correlations with extensive scans 
along (1/2, 1/2, L). As background we subtracted the in- 
tensity remaining at 80 K as there is no significant tem- 
perature dependence in the quasielastic scattering above 
80 K. The results presented here have been confirmed 
with polarized neutrons (see later) which is sensitive only 
to the magnetic cross-section. Typical scans through the 
correlated magnetic scattering are illustrated in Fig. [2l 
The momentum broadening is much larger than resolu- 
tion for the in-plane magnetic scattering at the (1/2, 1 /2, 
2) position and is similar at 2.5 and at 25 K. The strong 
elastic magnetic scattering observed in YBCOg.ss con- 
trasts with its absence in more highly doped Ortho-II 
YBCOe.s^ 



B. Correlation Lengths and Integrated Intensity 

To extract quantitative information we fitted the data 
to a Lorentzian squared profile convoluted with the four- 
dimensional instrumental resolution function. The fol- 
lowing profile gave a good description of the in-plane q- 
dependence of the elastic commensurate magnetic scat- 
tering as illustrated in Fig. O 
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Here A is an amplitude, Q = (27r/a)(ff, if ), Qo = 
(27r/a)(0.5, 0.5), and f is a correlation length that is 
isotropic within the a — b plane. As we shall see the 
out-of-plane correlations decay quite slowly along L so 
that modeling the in-plane correlations with only the 
two-dimensional vector, Q, is valid. As the sample is 
twinned, we use tetragonal coordinates for the remainder 
of the paper. Since a graphite analyzer was used only a 
limited integration over energy was done. However, since 
the central mode part of the spectrum is much narrower 
in energy at the lowest temperatures and stronger it fol- 
lows that ^ is to be interpreted not as the instantaneous 
correlation length but as one for near-static fluctuations. 



5 



This profile differs from the Lorentzian of our previous 
cold neutron work. In the present analysis the energy 
resolution was coarser (E/ = 14.6 meV with meV 
full width at half maximum) so that we could observe the 
scattering up to higher energies and temperatures where 
it is weaker. 

The Lorentzian squared profile follows from mean 
field theory for subcritical fluctuations.— The finite 
low-temperature correlation length can arise when ran- 
dom fields dominate over thermal fluctuations as pre- 
viously applied to the spin-glass region of LSCO and 
in two-dimensional magnets4^ Both a Lorentzian and a 
Lorentzian squared give a good description of the data 
but the correlation range extracted is naturally shorter 
for the Lorentzian squared. We note that in a previous 
work we used a Lorentzian squared to fit the correlations 
along the c-axis and a Lorentzian to fit the correlations 
within the a-b plane.'^S The Lorentzian squared has the 
advantage of having a finite two-dimensional integral over 
momentum unlike a single Lorentzian. 

From fits to scans along the [110] direction, convoluted 
with the resolution function, we extract an in-plane cor- 
relation length ^a6= 13 ± 2 A for the slow quasielastic 
fluctuations at 2.5 K. Fits to the data along the [001] 
direction gave using the coupled bilayer 

model developed previously, which also corresponds to 
a Lorentzian-squared profile^-i^ The in-plane and out-of- 
plane correlation lengths are the same within the lim- 
its of the accuracy. Moreover we find that the width in 
wave vector of the elastic scattering is largely indepen- 
dent of temperature up to at least ~ 50 K an example 
of which is shown in Fig. [2] for 25 K. The correlation 
lengths do not diverge with decreasing temperature as in 
insulating La2Cu04.^^ We have therefore the interesting 
result that the central mode spin correlations are nearly 
isotropic and temperature independent. The correlation 
length has a geometric rather than thermal origin, for it 
matches the hole spacing for a two-dimensional system, 
I = [a/^Jp) A~ 15 A where a is the lattice constant 
and p ^ 0.06 is the hole doping. The in-plane correla- 
tion length is shorter than the indirect estimate of the 
mean free path (average of / ~ 38 A for YBCOg.ss) from 
optical conductivity-- but comparable to our Lorentzian 
analysis.— Our results are similar to lightly doped and 
non-superconducting LSCO in its spin glass regime for 
concentrations of ~ 0.02-0.04 for which 1/^ also sat- 
urates at low-temperatures and increases only at very 
high temperatures 4i The difference here is that the short 
length scales coexist within the superconducting phase, 
and there is no observable incommensurate modulation. 
Incommensurate modulation has been observed in LSCO 
and YBCO and has been explained both in terms of 
stripes and quasiparticle excitations . 
Based on the half-width at half maximum of our elastic 
peaks, we place an upper limit of ~ 0.02 (defined as the 
peak splitting along [110]) for any possible incommensu- 
rability. We note that an incommensurate wave vector 
of order ~ 0.035 rlu (along [110]), as seen in supercon- 



ducting LSCO (Ref. [T^ for a similar hole doping would 
have been readily detected with the resolution of this ex- 
periment (see Fig. [2]). Though the incommensurability 
in YBCOe+a; may scale differently with Tc than in the 
monolayer LSCO system, our data could be described in 
terms of four peaks displaced from (tt, 7r)if broadened in 
momentum. 

Scans along the [001] direction with E^=14.6 meV are 
presented in Fig. [3]at 2.3 and 30 K in panels a) and b). A 
featureless background scan at 80 K has been subtracted. 
The solid line is a fit to Lorentzians centered around each 
integer L values such that an intensity could be obtained 
for each L. The measured integrated intensity at each L 
value is plotted as a function of the calculated intensity 
in panel c). The calculated intensity is based on the for- 
mula I{L) = /2(l/2,l/2,L)sin2(7r L(l - 2z)) which is 
the product of the anisotropic Cu^+ form factor and the 
bilayer structure factor with z — 0.36i^ The measured 
integrated intensity at each integar L value was corrected 
for resolution effects as outlined elsewhere4^ Good agree- 
ment is obtained except for the L=2 peak which has too 
low an intensity. While we do not currently understand 
why the L=2 peak is so weak, we note that on tighten- 
ing the energy and momentum resolution considerably 
by going to E/=5.0 meV, the ratio between the L=2 and 
L=l peaks agrees perfectly with theory as does the ratio 
between L=2 and L=3f2^ We attribute this anomalously 
low intensity to a resolution effect. Even though peaks 
appear at the integer positions as expected for bilayers 
coupled ferromagnetically from cell to cell, the peaks are 
broad at all temperatures. At high temperatures (30 K), 
the integer modulations weaken leaving two peaks cen- 
tered on Ij—1.7 and 5.1, the fingerprint of the bilayer 
structure factor. This indicates the effective absence of 
correlations between cells. The c-axis correlation length 
of 8 A may be thought of as a mean square transverse 
slippage by '--^ 2 A of the in-plane spin correlations from 
cell to cell along c. 

The diffraction pattern in Fig. [3] leads to important 
conclusions. First it shows that the spins in the bilayer 
sheets are antiferromagnetically coupled at all tempera- 
tures, i.e., the phase across the bilayer is (/) = tt. Any slid- 
ing of the spin density in one layer relative to the other 
would be clearly visible. Likewise any systemic phase 
difference between cells along c away from the in-phase 
= 0, would have caused diffraction peaks at non-integer 
L that are not seen. Thus the regions of correlated spin 
density in one cell maintain a ferromagnetic in-phase re- 
lation from cell to cell. The correlated spin region in 
one bilayer thus lies above the next cell along c. This 
is precisely the form of spatial fluctuation that precedes 
a critical quantum transition to 3D antiferromagnetism. 
The correlation range is short and strongly sub-critical 
but the pattern is highly organized and not random. If 
the spins were disordered in clusters as claimed — the 
regular integer-centered pattern would not be observed. 
Neither disordered localized clusters, nor a cluster spin 
glass as inferred from muon studies, are in accord with 
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FIG. 3: a) and b) The elastic magnetic correlations along the 
[001] direction at 2.3 K and 30 K with a background at 80 
K subtracted. The correlation length along [001] is ~ 1 unit 
cell. The 3D correlations weaken above 30 K and arise from 
independent bilayers. c) The measured integrated intensity 
(defined in the text) is plotted as a function of the calculated 
integrated intensity for each labelled L value. 



YBCO , Temperature Dependence 




150 200 250 300 



S~1200 
'c 



3 b) 80 K Subtracted Intensity 



X2 
< 



800 



• E=14.6 meV 
E =5 meV 



CO 

£1 
CD 

c 600- 



400 



200 



<2 



20 40 60 80 

Temperature (Kelvin) 



100 



experiment. The spins clearly smoothly reorganize into 
a pattern where the susceptibility at 3D AF momenta 
{L = n) is a maximum. The frozen 3D spin pattern we 
observe shows that a single phase exists in which super- 
conductivity and glassy spins coexist. 

At higher temperatures than those in Fig. [31 the 3D 
correlations become unobservable. The elastic scattering 
becomes part of the single inelastic response discussed 
later in which the scattering is from independent bilay- 
ers. Thus, both energy and temperature destroy the 3D 
correlations. 

The saturation of the correlation length at low- 
temperatures is very similar to that observed in lightly 
doped LSC04i The very short correlations both within 
the plane and along [001], indicate that YBCOg.ss is far 
from any critical point. Our results also differ from the 
expected behavior near a quantum critical point where 
the correlation length should continuously diverge, not 
saturate at a high temperature value. Even though fur- 
ther confirmation is needed at even lower oxygen doping, 
our results suggest that YBCOg.as is far from any quan- 
tum critical point. One possibility is that the transition 
out of the superconducting phase leads to a phase with- 
out long-range ordered antiferromagnetism. The possi- 
bility that the transition from superconductor to a Neel 
ordered antiferromagnet may be first order is unlikely 
because of the smooth and continuous metallic thermal 
conductivity across the critical dopingi^ 



FIG. 4: a) The temperature dependence of the peak magnetic 
intensity at huj=0 with E/=14.6 meV at Q=(l/2,l/2,2.0). 6) 
The temperature dependence of the elastic magnetic scatter- 
ing measured with E/=14.6 meV and 5 meV (C5 Chalk River) 
with a background measured at 80 K subtracted. The wave- 
vector lies at the maximum of the broad spin correlations at 
Q=(l/2, 1/2, 2.0) The different fixed final energies probe the 
magnetic fluctuations with different energy resolutions and 
hence different timescales. 



C. Temperature Dependence 

The temperature dependence of the elastic peak mag- 
netic intensity was studied in momentum scans along the 
[110] direction through Q=(l/2, 1/2, 2.0). The results 
would be similar for non-integer L as this is not a Bragg 
peak. As the elastic correlation length is known to be 
temperature independent (see previous section) the peak 
height of these scans represents the resolution integral 
over the spectral profile of the central mode at least as 
T— >0 where it is resolution limited. The results of this 
analysis are presented in Fig. HI The directly measured 
magnetic intensity as a function of temperature at tko=0 
and at Q=(l/2,l/2,2.0) is illustrated in panel a). The 
low temperature data was derived from scans along the 
[110] direction, examples of which are shown in Fig. [2l 
High temperature data, where the central mode is weak, 
were extracted from inelastic constant-Q scans by ex- 
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trapolating the magnetic intensity to hu!=0 as in Fig. 
[TT] Overlapping data taken at 80 K give good agree- 
ment between the two methods of extracting the elas- 
tic scattering in the absence of a strong central peak. 
Since the elastic magnetic intensity shows little change 
above 80 K, in panel 46) we plot the magnetic intensity 
at Q=(l/2,l/2,2.0) using 80 K as a background. For two 
different final energies, E/=14.6 meV and 5 meV with 
energy resolutions (full width at half maximum) of 1.5 
meV and 0.1 meV respectively, we find a decrease in the 
onset temperature on narrowing the energy resolution. 
This demonstrates that very slow spin correlations oc- 
cur as in metallic spin-glasses which show a systematic 
lowering of the onset temperature as energy resolution is 
tightened.— 

The gradual onset of magnetic scattering in Fig. [4] 
contrasts with that observed in the insulator where a 
sharp well-defined Neel temperature is observed;^ The 
growth of scattering with decreasing temperature might 
be thought to arise because more of the spin spectral 
weight, described by a single relaxation rate at each tem- 
perature, slows so as to fall within the energy resolution. 
This is an inconsistent view. If the central mode spectral 
profile were to correspond to a single Lorentzian relax- 
ation rate, j{T), that slows on cooling, then its inten- 
sity would fall to half its zero temperature limit when 
7(r) matched the resolution half-width. For the 14.6 
meV data this would indicate that 7(r) — 0.75 meV 
at T=20 K. Evidence against a single Lorentzian relax- 
ation rate, independent of the T-dependence of the static 
susceptibility, comes from cold neutron measurements^ 
where we know that the energy half width near 20 K is 
only 0.05 meV, not the much larger 0.75 meV width that 
would be required for a single relaxation rate model of 
the 14.6 meV data. As in a glass it is incorrect to as- 
sociate either the half- height intensity or the apparent, 
sensitivity-dependent, high temperature onset with any 
single relaxation rate. The growth of central mode inten- 
sity seems to correspond to the appearance with cooling 
of a hierarchy of longer spin timescales. 

Measurements at longer timescales and hence even 
finer energy resolution, ~ 5000 times finer, have been 
made with fiSH and have indicated a lower onset tem- 
perature of statically ordered moments i^i^ However the 
temperature where the peak has fallen to half its height 
has not decreased by such a large factor, but only to 
«10 K. The much slower decrease in onset temperature 
with decreasing energy resolution indicates that the elas- 
tic correlations are neither statically ordered, nor fluc- 
tuating with a single relaxation rate that slows with de- 
creasing temperature. 



D. Polarization Analysis 

We have verified the method of deriving the structure 
of the local moments from unpolarized data by measuring 
the polarized neutron cross-section since it is sensitive 
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FIG. 5: Spin-flip cross-sections for horizontal field (HP) and 
vertical (VF) geometries measured at huj=0. For all values 
of L, the magnetic scattering in the HF channel is twice that 
in the VF channel, indicating no preferred orientation of the 
moments thus corresponding to a paramagnetic phase. 



only to the magnetic cross-section. 

In Fig. [5] the solid line shows a fit constrained so that 
the HF is twice the VF intensity. The curve gives an ac- 
curate description of the data (the errors are about the 
size of the symbols). Since the HF spin flip intensity is al- 
ways twice the VF spin-flip intensity, it follows that there 
is no preferred orientation for the quasi-static moments. 
The in-plane spin correlations along the vertical [1, 1, 
0] z-direction were determined in two ways as shown in 
Fig. [51 The upper panel shows the HF-VF intensity in 
the spin-flip channel, which gives the S^z spin scattering. 
It would amount to half of the total magnetic intensity 
for a paramagnet. Because the nuclear contribution is 
normally large, the non-spin flip intensity is not usually 
used to determine magnetic cross-sections. Nonetheless 
in the lower panel we show the subtraction VF-HF in 
the non-spin flip (NSF) channel. The HF intensity con- 
tains no NSF magnetic scattering and is a measure of the 
(nuclear) background. The resulting non-spin-flip differ- 
ence, VF-HF, measures only one of the two components, 
the same Szz of the magnetic moment perpendicular to 
Q. It carries larger error because of the subtraction of 
two large nuclear signals. The excellent agreement and 
consistency we obtain between spin-flip and non spin-flip 
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FIG. 6; Subtraction analysis of the polarized scans along 
[001]. The upper and lower panels both measure quasielas- 
tic magnetic intensity along the [1, 1, 0] or z spin direction 
but determined by two different methods. The upper panel 
is from the spin-flip channel and gives Szz, since for y per- 
pendicular to z and to Q, HF^S'^^ 4- Syy and VF=5'yj, for all 
L. The lower panel is from the non spin-flip channel and, us- 
ing the HF non-spin-flip intensity as background, it also gives 
Szz, the fluctuations along the vertical z direction. 



methods adds credence to the analysis and experimental 
method. 

Scans along the [110] direction through the correlated 
scattering are presented in Fig. [7] for energy transfers 
of 0, 2.1, and 4.1 meV respectively. The fitted lines are 
constrained such that the intensity in the HF channel is 
twice that in the VF channel. The satisfactory descrip- 
tion of experiment illustrates that the scattering cross- 
section is paramagnetic at all energy transfers. There is 
no preferred direction for the spins, neither for the excita- 
tions nor for the static correlations. The magnetic cross- 
section can then be extracted at all energy transfers by 
subtracting the VF from the HF scattering. The param- 
agnetic cross-section at all energies is different than that 
observed in insulating YBCO where a significant out-of- 
plane anisotropy causes an energy gap of about 5 meV.^^ 

Our polarized neutron analysis of the elastic magnetic 
scattering proves that the spin orientation (polarization) 
of the static moments is isotropic. This is very differ- 
ent from the AF ordered insulator where the moments 
are pointing within the a — b plane. The spins in the 
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FIG. 7: Constant energy scans through the correlated mag- 
netic peak in the superconducting state at 9 K. The HF cross- 
section (solid circles) is seen to be twice the VF cross-section 
(open circles) at each energy transfer. The solid (broken) 
lines are Gaussian flts to the HF (VF) data that confirm that 
dynamic as well as static fluctuations are isotropic in spin 
direction. 



superconductor form a paramagnetic structure that we 
call a hedgehog phase. The isotropic polarization of the 
spins is not surprising given that the correlation lengths 
are nearly isotropic with the in-plane close to the out-of- 
plane correlation length. We then only expect significant 
anisotropy to develop when the correlation length grows 
to more than a few lattice constants in the a—b plane and 
along the c direction. Our results may be consistent with 
theoretical predictions of Hasselmann et al. (Ref. IstI ) 
and Lindgard (Ref. HI) based on spiral spin states. They 
predict the lightly doped region of the cuprate phase dia- 
gram is characterized by a significant dipolar frustration 
of the antiferromagnetic order leading to spiral correla- 
tions. The frustration is caused by holes (which reside 
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FIG. 8: The complete low-temperature inelastic and elastic 
spectrum obtained with polarized neutrons. To obtain the 
purely magnetic scattering we have subtracted the VF from 
the HF as discussed in the text. There are two spectral com- 
ponents, one centered around huj=0 and another associated 
with the broad inelastic scattering extending to higher ener- 
gies. 



primarily on the oxygen atom) causing Cu^'*' spins lo- 
cally to be ferromagnetically coupled. Earlier strong fer- 
romagnetic exchange through holes was found to cause 
chiral spin rotations and destroy order Models for the 
dispersion by Lindgard^S could result in an out of plane 
component to the moment as observed in YBCOg.as. 

Fig. [8] shows the complete low-energy spectrum of 
YBCOe.ss at 9 K using polarized neutrons. The magnetic 
cross-section for energy transfers greater than 2 meV was 
extracted from constant energy scans along the [110] di- 
rection over the ridge of scattering. The two channels 
were then subtracted and a Gaussian fitted to the HF- 
VF data. For for energies less than 2 meV a single point 
was measured at (1/2, 1/2, 2.0) in both the HF and VF 
channels and subtracted. The spectrum clearly shows 
the spins fluctuate on two time scales, one giving a sharp 
central elastic peak and the other a broader component 
peaked around 2 meV and smoothly decreasing with en- 
ergy up to ~ 35 meV (the highest energy measured). The 
spectral profile agrees with that derived from unpolarized 
neutrons as shown in Fig. [T] and as discussed in the next 
section. 
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FIG. 9: Out-of-plane and in-plane dynamic magnetic corre- 
lations showing the momentum dependence of the magnetic 
scattering at T=9 K for various energy transfers. The solid 
line for the scan along (1/2, 1/2, L) gives the behavior ex- 
pected for independent bilayers. The solid lines through the 
(H, H, 2) scans are results of fits to Gaussians. The horizon- 
tal bar represents the momentum resolution along the [110] 
direction. 



IV. INELASTIC RESPONSE 

In the previous section, we discussed the temperature 
and momentum dependence of the elastic and low en- 
ergy magnetic scattering. The lowest energy fluctuations 
within the central peak arise from spin correlations that, 
for the lowest temperature, become quasi-static on the 
neutron timescale. While their correlation lengths within 
the basal a-b plane and along the c-axis are short, these 
slow fluctuations are 3D in nature. In this section, we dis- 
cuss the inelastic scattering from faster spin fluctuations. 
We will show that the fast fluctuations are primarily two- 
dimensional in nature and will describe the temperature 
dependence of their profile in energy and momentum. 



A. Momentum profile 

The momentum dependence of the correlated magnetic 
scattering is illustrated in Fig. [9] with scans along the 
[001] and [110] directions. The solid lines through the 
scans along the (H, H, 2) direction are results of fits to 
a Gaussian of fixed width convolved with the resolution 
function. The good agreement shows that the in-plane 
correlation length is largely independent of energy just 
as earlier we showed its independence of temperature. 
The solid line through the scan along the (1/2, 1/2, L) 
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FIG. 10: a) The half width at half maximum in momen- 
tum as a function energy transfer at 4 K. The dashed line in 
the left-hand panel illustrates the spin-wave velocity (hence 
E=?i,C75=5907g) determined at high energy (Ref. |38|). b) — d) 
The temperature dependence of the half width at half max- 
imum at energy transfers of 2.1, 12.4, and 24.8 meV. The 
widths are from fits of a Gaussian profile to scans along (H,H) 
direction around Q=(l/2,l/2). 



direction is the result of a fit to the following equation 
representing the acoustic bilayer structure factor: 



S{L) = Af{L) sin2(7r L{1 - 2z)) + B. 



(2) 



Here A is an amplitude, (1 — 2z) represents the bilayer 
spacing with z=0.36, P{L) is the anisotropic Cu^+ form 
factor, and B is a constant. This model with two free 
parameters gives an excellent account of the data over a 
wide range in L. The spin excitations therefore arise from 
independent bilayers. Note that the optic fluctuations 
play no role as they lie above an onset of ~ 50 meV (Ref. 
1381 ). well outside the energy range of this experiment. It 
is interesting that the peaks which appeared at integral 
values of L for the quasielastic central peak are no longer 
observable. Therefore, the 3D behavior is only exhibited 
by the slow dynamics of the central peak, whose large 
susceptibility enhances 3D correlations, and not by the 
faster timescale of the inelastic channel. 



The observed in-plane magnetic scattering is commen- 
surate. There is no sign of the flat-topped structure nor 
of the hour-glass dynamic incommensurability observed 
at larger doping^ in Ortho-II YBCOg.s, nor of the well- 
defined incommensurate features of lightly doped normal 
and superconducting LSCO^ and LBCOi^ Instead Fig. 
|9] shows a simple in-plane momentum broadening that is 
very weakly dependent on energy between 2.1-24.8 meV 
with an average Gaussian-fitted half width at half max- 
imum (hwhm) (Fig. ^ of 0.09 ± 0.02 A~\ Indeed the 
scans at various energies in Fig. [9] are well described by 
a constant width. The very small increase of momentum 

broadening with energy (0.0.09 at 3.7 meV to 0.10 A 
at 16.5 meV) mirrors the slope from the large spin-wave 
velocity which is known from high-energy measurements 
on the same sample^® and is indicated by the dashed 
line in Fig. [TO] a). In contrast, for larger doping the mo- 
mentum broadening decreases with energy because of the 
hour-glass dispersionj^ For E=0 the Gaussian half width 

(hwhh) is 0.06 ± 0.01 A or about two thirds of the 
broadening of the spin excitations. 

Apart from momentum broadening, the commensurate 
profile along [110] resembles that in the insulating region 
of YBCO and LSCO despite the fact that YBCOg.ss is 
metallic and superconducting. For systems with a highly 
suppressed Tc such as YBCOg.as a dynamic stripe model 
may therefore not apply. It is possible that, for very low 
concentrations of holes, regions of correlated spins are fa- 
vored over a quasi-one dimensional alternation of spins 
and charge. We note, however, that the expected value 
of the incommensurability (denoted as 5) is small. Based 
on the monolayer LSCO systeni^, where the incommen- 
surability for 6 % doping is ^ 0.05 r.l.u. or 0.08 A 
along [10 0], we would expect for our sample an incom- 
mensurate wave vector of 0.06 A ^ or 0.025 in H as our 
scan along [H H 0] passes between putative incommensu- 
rate peaks. This is close to the hwhh of 0.06 ± 0.01 A 
or 0.025 in H we measure for E=0 correlations (Fig[2]). 
Since the elastic fwhm resolution in H is 0.01, an incom- 
mensurate wave vector of 0.025 would have been easily 
resolvable if it were well-defined. Thus the absence of 
incommensurable peaks stems not from the more than 
adequate resolution but from sources of intrinsic broad- 
ening such as the short range of the oxygen chain order 
or the dominance of magnetism over superconductivity 
when Tc is suppressed. 



B. Energy spectrum 

We determined the energy dependence of the magnetic 
scattering from constant Q and from constant energy 
scans for energy transfers up to ^ 25 meV. The peak 
magnetic intensity was extracted from the constant Q 
scans on the ridge at (II,II)=(0.5, 0.5) from which a back- 
ground averaged between (0.3, 0.3) and (0.7, 0.7) was 
subtracted. An independent determination of the mag- 
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FIG. 11: The magnetic scattering as a function of energy 
transfer is plotted for 4, 80, and 295 K. The data was 
obtained using unpolarized thermal neutrons with a non- 
magnetic background subtracted. The sohd line is the result 
of a fit to the modified Lorentzian lineshape described in the 
text. The data around hto =0 has been removed so that the 
purely inelastic spectrum can be analyzed free from the cen- 
tral quasielastic component. The horizontal bar is the energy 
resolution at the elastic line. The elastic scattering around 
hLj=0 has been removed. 



netic signal was obtained at selected energies in scans 
over the ridge along [110] direction such as those in Fig. 
ini The peak magnetic intensity from a Gaussian fit gave 
excellent agreement with the constant Q method. 

The magnetic continuum spectrum thus determined 
exhibits spin relaxation (Fig. [TT|l that broadens and 
weakens as temperature rises. At low temperature the 
scattering is peaked at around 2 meV, the limiting relax- 
ation rate for the spin excitations for the short-range or 
glassy spins. The 2 meV relaxation rate lies much below 
any energy derived from the product of momentum width 
and high spin velocity. At 295 K the scattering is broad 
and almost energy independent. The absence of data 
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FIG. 12: The temperature dependence of the linewidth 2r 
and the susceptibility are plotted. The broken line in panel 
is the result of a fit to a straight line through the origin. 
The broken line in panel b) is the result of a fit of the high 
temperature data to l/T. This temperature behavior breaks 
down for temperature below ~ 50 K. 



around ~ 20 meV is caused by a strong phonon which 
is difficult to correct for with unpolarized neutrons. We 
know, however, that the magnetic scattering continues 
smoothly through the phonon region from the polarized 
neutron analysis (Fig. [8]). The spectrum differs consid- 
erably from that observed in YBCOe.s (Tc=59 K) at low 
temperatures (Fig. [1] which rises to a resonant peak in 
the susceptibility at 33 meV-. By suppressing the critical 
temperature from 59 K to 18 K, the spectral weight has 
shifted to low-energy at temperatures well below the on- 
set of superconductivity. Instead of the positive slope of 
the dynamic susceptibility for the 59 K superconductor as 
the energy rises to a 33 meV maximum at a resonance, 
we find a negative slope corresponding to a weakening 
susceptibility as the energy rises beyond the maximum 
of 2 meV. These two very different behaviors require a 
major electronic reorganization to have taken place. 

To obtain the relaxation rate and amplitude of the 
spectrum at Q=(l/2,l/2) and below ~ 25 meV we fit- 
ted the scattering, with the quasielastic central peak re- 
moved, to the following modified Lorentzian, which cor- 
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responds to overdamped excitations with a single energy 
scale, the relaxation rate F, 

5(Q,u;) = x'Wc^) + l]p2^- (3) 

where x' is the static susceptibility, and [n{ijj) + 1] is the 
Bose factor.. The solid lines in Fig. [TT] shows that this 
model with only two adjustable parameters at each tem- 
perature gives an excellent description of experiment over 
a very broad range in temperature and energy. The fast 
spin dynamics therefore relax with a single temperature 
dependent time constant. The amplitude and relaxation 
rate are shown in Fig ll2l up to room temperature. 

As illustrated in Fig. [1] the spectrum of spin fluctua- 
tions at low temperatures in YBCOg.as (Tc=18 K) differs 
considerably from that observed in Ortho-II YBCOg.s- 
YBCOe.s displays a clear resonance peak at ~ 33 meV in 
the superconducting phase. High-energy measurements 
in YBCOg.s have shown the resonance defines a cross- 
over point between low-energy anisotropic incommensu- 
rate fluctuations and high-energy spin-wave-like fluctua- 
tionsi^ In contrast to this hour-glass shaped dispersion, 
in YBC06.35, no clear resonance peak is observed. Even 
though the neutron scattering response peaks at ~ 2 meV 
the spin response is relaxational, and does not scale with 
Tc as does the sharp resonance peak in YBCOe.s. The 
relaxational feature does carry, however, a substantial 
fraction of the low-energy spectral weight. 

The absence of a spectrally sharp resonance can be 
attributed to several effects. One possibility is that spin 
disorder introduced through short oxygen chains strongly 
damps the resonance peak. We find this conclusion un- 
likely, for if the chain disorder were to introduce such a 
large effect as to completely destroy a resonance peak, 
we would expect a strong effect on T^ (which is mea- 
sured to be very well defined in our crystals). We would 
also expect that, in a completely disordered sample, the 
resonance peak would be absent. This is clearly not the 
case as a resonance peak is observed in oxygen disor- 
dered samples near the YBCOg.s concentration and also 
in heavily overdoped materials. ^'^ While it has been found 
that impurity doping affects the width of the resonance 
peak, the ratio of the resonance energy to that of Tc 
remains unchanged with impurity dopingj^ Scaling the 
resonance energy with Tc, from our YBCOg.s results^ we 
would expect a resonance peak at 10 meV in YBCOg.ss, 
clearly not the case in our measurements. Therefore, 
chain disorder cannot completely account for the absence 
of a strong resonance peak in YBCOg.as. Spectral broad- 
ening is known to occur, for example in YBCOg.s, where 
disorder increases the resonance width at low tempera- 
tures from 10 meV in ordered Ortho-II to 25 meV in a 
disordered sampleJ^ It has been found that while impu- 
rity doping increases the width of the resonance peak, 
the ratio of the resonance energy to that of Tc remains 
unchanged with impurity dopingi^ The 2 meV feature 
may possibly be viewed in terms of depression by large 



damping near the critical doping as the spin response 
evolves into critical scattering. 

A low resonance energy might also arise because the 
incommensurate wave vector is small or absent. If we 
adopt the stripe model^'^, where the resonance is a cross- 
over energy that separates incommensurate fluctuations 
from spin- waves, we would then expect a very low energy 
resonance because of the low incommensurability. 

Our observation of strong but only relaxational spin 
excitations in YBCOg.ss indicates that any resonance 
does not need to be well-defined for superconductivity to 
occur. However spectral weight from the resonance re- 
gion can be important in heavily doped systems such as 
YBCOg.gs where it is the only easily detectable magnetic 
spectral feature and where it gives more than enough 
change in magnetic energy to account for the supercon- 
ducting condensation energy.— The belief that a reso- 
nance cannot account for superconducting pairing exists 
because the resonance comprises so little of the total 
spectral weight (Ref. f65h . Other theories link the res- 
onance to the pairing boson (Refs. [sl, 113) ■ The de- 
coupling of the resonance energy from Tc is illustrated 
by LBCO at a doping of 1/8 where the resonance en- 
ergy occurs at ~ 60 meV despite the fact that Tc is al- 
most zerojii Thus the resonance energy may scale with 
Tc in YBCO over a limited range of hole doping, but it 
does not near 1/8 doping nor for different cuprates.^^ In 
LSCO two peaks have been observed in x" ^^^d the lower 
energy peak is not obviously connected with cross over 
from stripe-like fluctuations to spin waves.®* Instead the 
resonance energy may scale with the incommensurability 
that is driven by accommodation of holes rather than by 
pairs, and hole doping is small for YBCOg.ss. 

The susceptibility x! and full-width 2F as a function of 
temperature are plotted in Fig. [121 Panel h) shows that 
the susceptibility on cooling initially grows strongly as 
1/T (broken line), then slows below 50 K and finally sat- 
urates below ~ 10 K. Panel a) shows that the linewidth 
(full width at half maximum) decreases linearly with tem- 
perature until it saturates below ~ 50 K. The broken line 
in Panel is a high temperature fit to 

2r= (1.40±0.15)fcBr, (4) 

where F and /csT are given in meV. The proportionality 
of the relaxation rate to temperature implies that tem- 
perature is the dominant energy scale above ~ 50 K. As 
shown in the next section, uj/T scaling is expected and 
found to be obeyed over a broad range in temperature 
and energy. 

For critical scattering in a two dimensional periodic 
Heisenberg system, the energy linewidth (F) is related to 
the correlation length {£) by KT — Tic/£^ where he is the 
spin-wave velocityi ®^i^° Saturation of the characteristic 
energy linewidth would then be attributed to saturation 
of the correlation length. Instead, as discussed above, 
the spatial correlation length determined from constant- 
energy scans is almost independent of temperature below 
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50 K (Figs. mUSl). This indicates that the length is ge- 
ometrically limited, but its magnitude is not linked by 
velocity to the dynamic width. For example the 4 K cor- 
relation length of l/,f = 0.09 A times the velocity of 
590 meV-A gives an energy of 50 meV, much too large 
to account for the 2 meV relaxation rate. Thus any re- 
lation between spatial and temporal scales breaks down 
for the spin relaxations and more so for the spectrally 
narrower central peak. Instead measurement shows that 
r saturates at the same temperature where the spectral 
weight transfers to the much slower central peak, as dis- 
cussed later, with an almost constant correlation length 
for both. A similar result was obtained in La2-a:Sr3:Cu04 
where the correlation length saturates at low temper- 
atures for a similar hole doping. ''^ Any effective veloc- 
ity linking lengths to times would have to be very small 
as might occur for the near localized excitations of non- 
periodic glassy spins. From the temperature dependence 
of its dynamics, YBCO6.35 cannot be interpreted as being 
close to a quantum critical point. 

The saturation of the spatial and dynamic scales re- 
lates to a recent theoretical t-J modcl^^ that predicts 
that F should decrease to zero upon entering the insu- 
lating antiferromagnetic phase. Whether F drops to zero 
continuously (in a second order manner) or discontinu- 
ously (if the quantum transition is first-order) remains to 
be investigated at lower doping. 

Perhaps what is most surprising in Fig. [12] is that no 
clear effect of the superconductivity is observed in the ex- 
citation spectrum of YBCOe.as on passing through the 
superconducting transition. Neither does the elastic scat- 
tering change near Tc, as noted in Section III B. Both 
the inelastic and elastic response emerge continuously as 
the spins slow down with cooling. In contrast, in Ortho- 
II YBCOg.s a sharp and distinct signature of coherent 
superconducting pairing was observed causing a suppres- 
sion of the low-energy incommensurate intensity below 
huj 16 meV. Simply scaling this superconducting gap 
by Tc, we would have expected to see a suppression of the 
magnetic scattering below energy transfers of ~ 5 meV, 
clearly not the case in YBCOe.as based on our thermal 
and also cold neutron worki^ As pointed out later with 
regards to the total moment sum rule, the absence of spin 
suppression cannot be attributed to macroscopic phase 
separation for we observe the expected total amount of 
magnetic scattering. 

Our results may be associated with a microscopic de- 
coupling of magnetic and superconducting order param- 
eters. The a phase transition to a state where local 
magnetic fields coexists with superconductivity has been 
inferred by Miller et al. using /iSR by implanting the 
positively charged muon;^^ In the presence of the large 
carrier depression caused by the strong Coulomb screen- 
ing on a scale of ^1 eV, in what is a very low density 
system near a transition to localization, we find it likely 
that the field, and indeed the phase, seen by or caused 
by the muon may differ from that detected by more gen- 
tle non-Coulomb probes such as neutron scattering and 



NMR. 

One possibility is that locally antiferromagnetic re- 
gions are formed (within our measured correlation range) 
separated by metallic regions. Although different topo- 
logically, this has the same underlying antiphase-domain 
physics that causes the well defined stripes observed in 
the nickelates.'''^ A similar island structure has been pos- 
tulated for LSCO in the insulating spin-glass region. It 
could be the case that at higher hole dopings, a one- 
dimensional stripe structure is favored over islands which 
exist near the critical doping. A related model of stripes 
and superconductivity addresses the different kinds of 
hole states >2i 

Nonetheless our observation of a continuous evolution 
to a subcritical 3D spin pattern, without a transition in 
temperature, is at odds with a picture of physically sep- 
arated islands, or of a cluster spin glass, and certainly 
with phase separation, for which the highly organized 
3D spin pattern along [0 1] would not occur. In par- 
ticular the continuous transfer of spectral weight from 
the relaxational to the central mode, as discussed later 
(Fig. [TB)) . is a strong argument in favor of a single phase 
in which slow short-range glassy order and spin excita- 
tions coexist with superconductivity. The absence of a 
spin anomaly at Tc at low doping indicates that it is the 
holes rather than the coherence of pairs that break up 
the underlying antiferromagnetic spin pattern. Only at 
larger doping Tc = 59 K with stronger superconducting 
order, is the onset of pair coherence large enough to cause 
a superconducting gap in the spin response. 



C. uj/T Scaling 

In early studies x" in disordered YBCOg.s was found 
to follow a simple cj/T scaling relation over a broad 
range in energy transfers and temperaturesi^ The same 
scaling relation successfully describes the temperature 
dependence of the local susceptibility in underdoped 
La2-2:Sra,Cu04.^^'^^ The scaling relation, which assumes 
that the dominant energy scale is set by the tempera- 
ture, predicts the following temperature and energy de- 
pendence for the susceptibility, 

x-(c.,r^o) ^^"^"^"H~''^^-J- 

Here the term x"(a;,T = 0) represents the limit of x" 
as the temperature goes to zero. In Ortho-II ordered 
YBCOg.s the temperature dependence at all energies up 
to ~ 20 meV is well described by only the first term in the 
expansion with coefficient (ai). A similar expansion de- 
scribes the the excitations up to ~ 30 meV in disordered 
YBCOe.s.— Authors of recent studies of hole doped and 
electron doped cuprates have chosen to use another form 
of scaling analysis where x"iQAF^^, T) is plotted against 
Lo/T to find a single universal curvej^i^l This particular 
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FIG. 13: Constant energy scans for luo^lX meV and 24.8 
meV for a series of temperatures. The inelastic scattering 
is commensurate with no observable change in width as a 
function of temperature and energy transfer. The solid lines 
are the results of fits to Gaussians with a sloping background. 
The u/T scaling analysis was conducted based on Gaussian 
fits such as those represented here. 



analysis is a low-energy approximation and has been pre- 
sented in connection with the possibility of the close prox- 
imity of a quantum critical point. We have chosen to use 
the well established analysis which assumes only the less 
restrictive hypothesis that temperature is the dominant 
energy scale and so should apply to all energy transfers. 

To obtain an amplitude (and hence x") as a function of 
energy and temperature, constant energy scans were con- 
ducted along the [110] direction for a fixed temperature 
and energy transfer. The value of L was chosen based 
on the bilayer structure factor to maximize the scattered 
intensity. For our scaling analysis, temperatures in the 
range of 2 - 293 K were investigated and energy transfers 
in the range of 2.1 - 24.8 meV measured. Elastic scat- 
tering does not figure in the scaling analysis. Typical 
scans over a broad range in temperature and energy are 
presented in Fig. [T31 The correlated scattering is well 
described by the fitted Gaussians. We have also made 
use of the fact that the momentum width is found to be 
independent of energy and temperature. Each point in 
the scaling analysis is then based on a single amplitude 
parameter from the constant energy scan and is propor- 
tional to the local susceptibility. Recent results show that 
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FIG. 14: The normalized susceptibility x"(w, T)/x"(a;, 0) is 
plotted as a function of temperature at 2.1, 12.4, and 24.8 
meV. The solid lines are results of fits to the LojT scaling 
analysis described in the text. The fit describes the data over 
a broad range in temperature and energy. 



the momentum profile broadens and evolves at high en- 
ergies into two peaks whose velocity is the same as that 
observed in the cuprate insulator— but this occurs out- 
side the range hio < 35 meV of the current experiment. 

The normalized susceptibility (defined as 
x"{uJ,T)/x"{uJ,0)) is plotted in Fig. [H for three 
typical energy transfers. We then simultaneously fit the 
scaling equation of Eqn. [5] to all data as shown by the 
broken line in Fig. 1151 The entire susceptibility was best 
described by the inclusion of two terms in the scaling 
analysis: 
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FIG. 15: A plot of the normalized x" derived from constant 
energy scans for all energy transfers studied {hui — 2.1, 3.7, 
8.3, 12.4, 16.5, and 24.8 meV). The solid lines are resuhs of 
fits to the scaling equation discussed in the text. The a;-axis 
is plotted on a logarithmic scale. 



x"(^,r = o) 

^ arctan ^0.90 ± 0.04 + 2.8 ± 0.3 . (6) 

The inclusion of higher terms in the expansion did not 
noticeably improve the fits. The sohd hnes in Fig. [T^are 
typical predictions of the scaling equation. It is seen that 
u/T scaling gives an excellent description of the data over 
several decades. The value of ai is very close to the value 
of 1.1 obtained in the YBCOg.s Ortho-11 superconductor 
for temperatures above Tc and energies well below the 
resonance energy. The values for ai and a,2 differ con- 
siderably from those obtained in LSCO for similar hole 
doping^i While this implies a common temperature de- 
pendence to the magnetic fluctuations over a wide range 
in the underdoped YBCOe+x phase diagram, the scaling 
parameters do not necessarily universally apply to other 
cuprate superconductors. 

Even though scaling is approximately obeyed over a 
broad energy and temperature range, the lowest panel 
of Fig. [HI illustrates that the susceptibility at 2.1 meV 
does show evidence for a subtle breakdown below 50 K. 
We attribute this breakdown to the low-temperature re- 
gion where the energy full-width 2r ceases to decline lin- 
early with temperature. The reduction below scaling in 
the susceptibility at 2.1 meV below 50 K occurs at the 
temperature where the weight in the relaxational spin 



excitation begins to transfer to the quasielastic central 
mode as will be seen (Fig. [TB)) . 

Previously in lightly doped insulating LSCO, a simi- 
lar breakdown oi uj/T scaling was observed. The energy 
scale required to account for the breakdown was asso- 
ciated with an out-of-plane anisotropy which causes a 
gap in the excitation spectrum and hence a new energy 
scalei^ In YBCOg.ss we cannot attribute the breakdown 
of scaling to the presence of such an anisotropy as we 
observe the elastic correlations to be paramagnetic or 
isotropic, in contrast to the behavior of the insulator. 
However, the new energy scale required to account for 
the breakdown of scaling can be attributed to the con- 
finement of the spins as discussed in the context of cor- 
related regions in the previous section. 

In Ortho-II YBCOg.s scaling is obeyed in the normal 
state above 59 K and for energy transfers below ~ 20 
meV. The breakdown of scaling at high-energy transfers 
in ortho-II was attributed to the presence of a resonance 
peak at 33 meV introducing a new energy scale. The 
breakdown for temperatures below the onset to super- 
conductivity was characteristic of the superconducting 
gap opening and suppressing magnetic fluctuations be- 
low 16 meV. Neither a superconducting gap, nor a strong 
resonance peak is observed in YBCOg.as, and therefore 
scaling is observed over a broad range in temperature and 
energy. Theoretical motivation for scaling can be found 
in theories of marginal Fermi liquids (Ref. PtST) . close 
proximity of a quantum critical point (Ref. l2ll ). stripe 
liquid (Ref. 79), damping of a collective mode (Ref. 80) . 
and from a memory- function spin- wave theory (Ref. ISlI ). 

The fact that we observe scaling over a broad energy 
and temperature range in YBCOg.as shows that any sub- 
stantial spin pseudogap is absent. In contrast a charge 
pseudogap has been investigated for more highly doped 
YBCOg+a; by optical and NMR measurements which re- 
veal a suppression of 1 /TiT and of the optical conductiv- 
ity below 30 meV.*^''*'^ The presence of a spin pseudogap 
would imply that spin response x" would be suppressed 
at low temperatures if it were to follow the NMR relax- 
ation l/TiT or the optical conductivity. On the con- 
trary the staggered spin susceptibility and instead grows 
monotonically on cooling and we observe scaling and mo- 
ment conservation over a broad range in temperatures. 
Sutherland et ai recently suggested from transport mea- 
surements that the pseudogap energy scale may lie at ~ 
160 meV in low-doped materials. Therefore, the effects 
of the pseudogap are likely only felt at very high ener- 
gies as recently revealed by the large spin spectral weight 
loss above 160 meV in the high-energy magnetic excita- 
tions 



V. INTEGRATED INTENSITY AND SUM RULE 

To obtain a quantitative measure of the spectral 
weight, and of the fraction of the total moment which 
resides in the low-energy central peak, we consider the 
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FIG. 16: The integrated intensity in absolute units as a func- 
tion of temperature for energy transfers below 30 meV. The 
upper panel shows the transfer on cooling of weight from the 
broad inelastic relaxational Lorentzian to the central peak 
around hij=0. The lower panel depicts the total integral and 
shows that the sum rule is conserved at low temperatures. 
The decrease in total spectral weight at high temperatures is 
due to the broad energy width of the modified Lorentzian. 



sum rule for two uncoupled Cu02 layers. Ignoring the 
effect of the chains, and therefore taking only two Cu^+ 
ions per formula unit, the total spectral weight integrated 
over all energy and momentum is given by the total mo- 
ment sum rule: 



y y (fq S{q,uj) = 2 



1)9' 



(7) 



where J d'^q is the momentum space integral over the 
correlated peak. In this equation the factor of 2 comes 
from the fact that we are taking two Cu^+ ions per for- 
mula unit. Here we have assumed that S'(q, has been 
corrected for the bilayer structure factor and the Cu^"*" 
anisotropic form factor. Doubtless the total moment sum 
rule will be somewhat reduced as some of the spins will 



be destroyed by the doped holes. Since the doping is so 
small, p ~ 0.06 for YBCOg.as, and since the oxygen holes 
that screen the Cu^^ are not seen at neutron momenta, 
we expect the total moment sum rule to be approximately 
obeyed. Equating the above integral to the cross-section 
for paramagnetic scattering 



dcj / d3g[nM-f l]x"(q,c^) 



-fiWsis + i). 



(8) 



gives for 8=^ and g = 2 a, total integral of 1= 2/j,^. This 
limit provides a useful benchmark with which we can 
compare our integrated intensities. We determined the 
absolute intensity with two methods. All elastic and in- 
elastic intensities were calibrated to the known constant- 
Q integral over a transverse acoustic phonon at (0.15, 
0.15, 6) as described in Ref. l6l. This gives one value 
for the calibration constant. As an independent check, 
we also calibrated just the elastic scattering by measur- 
ing a series of Bragg intensities corrected for extinction, 
and fitting the calibration constant as described in the 
Appendix. Both methods agree well for the elastic scat- 
tering justifying the quantitative analysis used here. 

In YBC06.35, at 2 K, we observe a total moment below 
30 meV of 0.031 ± 0.006 The apportionment of the 
integrated spectral weight to the central peak and to the 
modified Lorentzian is given in the upper panel of Fig. 
[TBI The lower panel shows that the total integral (sum 
of central peak and modified Lorentzian) is conserved at 
all low temperatures below 150 K thus confirming the 
total moment sum rule (response above 30 meV~ 350K 
is not expected to change at these temperatures). The 
small reduction at high temperatures is a result of the 
energy scale of the modified Lorentzian becoming larger 
than the integration range of to 30 meV; it does not 
represent a loss of integrated spectral weight. 

As the temperature falls below ^ 50 K Fig. [11] shows 
that spectral weight transfers from the broad inelastic 
scattering to the central peak its gain being fully com- 
pensated by the loss of weight from the excitations. The 
conservation of spin weight is arguably the clearest ev- 
idence that we are dealing with a single phase system. 
The slowing of the relaxational response couples to and 
drives up the central peak as demonstrated by the excel- 
lent description provided by our coupled soft mode cen- 
tral mode modeli^. The model corresponds to inelastic 
excitations and a central peak that arise from the same 
single bulk phase. The extremely slow dynamics of the 
central mode can be viewed as the slow tumbling of the 
many correlated spins, and the fast inelastic relaxation 
as the excitations of these glassy spins. The central peak 
cannot be attributed to an impurity or macroscopic phase 
separation. Macroscopic phase separation has been cited 
as a possible reason for the coexistence of static magnetic 
order and superconductivity in La2Cu04+j,'^^, a view not 
supported by our data on high-quality crystals. 
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FIG. 17: The low temperature magnetic moment observed 
with elastic thermal neutron scattering as a function of the 
hole doping for both YBCO and LSCOi^ systems. The YBCO 
point at p=0.09 is an upper bound as described in the text.«4S 
The uniform suppression of the elastic magnetic scattering 
with no anomaly on entry to the superconducting state at 
Pc=0.055 indicates that holes not pairs destroy the moment. 



The transfer of spectral weight from the inelastic chan- 
nel to the central peak occurs at the same temperature 
as the saturation in the dynamic linewidth where we ob- 
serve a breakdown of cj/T scaling for low energies. The 
breakdown arises from the transfer to the gradually slow- 
ing down lowest-energy spins, forming a new energy scale 
represented by the central peak and not included in the 
scaling analysis. This conclusion differs from the scal- 
ing analysis performed on Li doped La2Cu04 which at- 
tributes the breakdown oi uj/T scaling to the close prox- 
imity of a quantum critical point. ^"'^ Nonetheless there is a 
critical point to antiferromagnetism in the vicinity and, 
despite very short range correlations, the central mode 
we see foretells its arrival at lower doping. 

In the YBCOe+a; system, the in-plane exchange con- 
stant is estimated to be J '--^ 125 meV implying that the 
spin excitations would extend to about Tilj ^ 250 meV. 
If the spectral weight were evenly distributed over all en- 
ergies we would expect to observe (30 meV/250 meV) x 
2 ~ 0.24 which is far above what is seen experimen- 
tally in YBCOe.as- Even recognizing that the density of 
states is weighted to high energies in a 2D antiferromag- 
net, we calculate that 0.06 of the spin spectral weight lies 
below 30 meV, a calculation known to be too large be- 
cause of the errors in spin-wave theory. We would then 
expect to observe 0.06 x 2 /i^ ~ 0.12 still much 



larger than the 0.03 we observe. To obtain a better es- 
timate of the expected total spectral weight, we rely on 
the spectral density of 4 ^\/eV obtained in insulating 
YBCOg.is'^ Taking the low-energy limit where the spin 
excitations disperse linearly with energy, the total ex- 
pected integrated moment should be '--^ 4 ^\/eV x 30 
meV /tt ~ 0.04 This agrees within error with the to- 
tal integral of 0.03 measured in YBCO6.35. The fact 
that we get good agreement with the total moment of the 
antiferromagnetic insulator provides additional evidence 
that the central peak is a necessary macroscopic property 
of single-phase YBCOe.as- 

The total integrated spectral weight up to 30 meV in 
YBCO6.35 agrees very well with that observed in more 
highly doped Ortho-II YBCOg.s despite the fact that the 
latter lies further from the antiferromagnetic insulating 
phase. As the integral from to 40 meV in ortho-II 
gave 0.05 /i^, we would expect that the integral from 
to 30 meV in YBCOg.ss should give ~ 0.04 /i^ in good 
agreement with what we measure. 

An important conclusion is that, within the YBCOg+x 
phase diagram, our results demonstrate that the total 
moment is conserved not only with temperature, but also 
with doping from the insulating region to deep within the 
superconducting phase. Therefore, the spectral weight 
is not removed on entry into the superconducting phase, 
but simply redistributed in energy and momentum. From 
YBCO6.35 to YBCOe.s, more spectral weight is moved 
from the elastic and 2 meV region to higher energies 
around the resonance energy at 33 meV (as illustrated in 
Fig. [T|). The apparent lack of spectral weight at larger 
doping may be associated with lack of visibility caused 
by a broadening in momentum and not necessarily by a 
removal of spectral weight 4^ A detailed analysis of the 
total moment and various contributions to the integral 
has been conducted by Lorenzana et al.^'^ 

In Fig. [17] we plot the dependence of the low tem- 
perature statically ordered moment on hole doping for 
both YBCO (see Ref. and this work) and LSCQiI. 
The moment includes the central peak, previously shown 
to be quasielastio^, and integrated over its momentum 
breadth. The ordered moment was obtained by calibrat- 
ing against Bragg peaks as described in the Appendix. 
For YBCO6.35 and YBCOe.s the ordered moment was ob- 
tained from a comparison with the nuclear Bragg peaks. 
The datum for YBCOe.s, where no magnetic Bragg peak 
is observable^^, is an upper bound placed using neutrons 
and has been confirmed by NMR.®^ With increase of dop- 
ing there is no break in the uniform suppression nor an 
anomaly as the superconducting phase is entered. We 
infer that it is the accommodation of holes by spin reor- 
ganization rather than pairing of charges that destroys 
the static moment. 

Significantly Fig. [T7] shows that the g-integrated 
quasielastic moment in YBCOg.as maps directly onto 
the curve for the measured incommensurate peaks in 
LSCO. It is surprising then that transport measurements 
at milli-Kelvin temperatures^* imply that LSCO and 
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YBCO have very different carrier densities for similar 
hole doping (LSCO is more localized). In contrast the 
static ordered moment on the THz timescale) indicates 
that LSCO and YBCO are very similar. This conclusion 
needs to be verified in YBCO to determine what differ- 
ences may exist between the LSCO and YBCO phase 
diagrams nearer the critical hole doping of pc=0.055. 



VI. DISCUSSION 

We have shown that the spin dynamics are governed 
by two distinct energy scales, one described by an energy 
resolution limited central peak indicating slow fluctua- 
tions which become static at the lowest temperature on 
the thermal neutron time scale, and a second fast com- 
ponent with a low temperature characteristic energy of 
2r=5 meV. The width of the fast component scales with 
temperature above 50 K and saturates at low temper- 
atures. We observe no well-defined resonant feature as 
occurs for larger hole doping. The momentum depen- 
dence illustrates that the correlations of the fast inelastic 
component are highly 2D in nature. 

The slow fluctuations, in contrast, are not described by 
a single relaxation rate and display 3D correlations at low 
temperature. The correlation lengths (determined from 
a Lorentzian squared) are short-ranged with S,ab— 13 ± 2 
A and $c=8 ± 2 AThe spins exhibit isotropic orientation. 
The quasielastic scattering corresponds not to static Neel 
order, but to a gradual slowing down of the dynamics 
with temperature. The slow component only becomes 
appreciable at low temperatures below fesT ~ 50 K ~ 
5 meV, where the dynamic width of the fast component 
saturates. It is not clear what sets the energy scale of 
the fast component but we regard it as the excitation 
spectrum of a short-range spin glass structure. 

The results cannot be attributed to the presence of 
macroscopic phase separation of a superconducting re- 
gion and antiferromagnetic regions with a range of hole 
doping. Our samples show a single sharp superconduct- 
ing transition which rules out the presence of other super- 
conducting domains with different hole doping. If phases 
with a range of Neel temperatures were present we would 
have observed resolution limited peaks at the antiferro- 
magnetic wave vector. Instead we observe an elastic cross 
section that is broad in momentum along the c-axis and 
also within the a — b plane. The most telling evidence 
for single-phase behavior is that the central peak gathers 
spectral weight as the inelastic channel loses weight. The 
two features are therefore dynamically coupled while con- 
serving total moment confirming an earlier hypothesisi^ 
This conclusion is further verified by the single lattice 
constant measured by X-rays in our samples. As shown 
by Liang et ai, the lattice constants (and in particular 
the c-axis) are very sensitive to hole doping in the heavily 
underdoped region of the phase diagram. "^^ We do not ob- 
serve a spread in lattice constants and therefore conclude 
that our sample is homogeneous. 



The work presented here is very different from early 
studies conducted in the underdoped region of the phase 
diagram. Work by J. Rossat-Mignod et al. did not in- 
vestigate the region of the phase diagram with a Tc less 
than 50 K and did not observe directly the two dynami- 
cally coupled timescales differentiated in the present mea- 
surements.*^ Their samples were not superconducting for 
oxygen content less than 6.40. More recent work by Mook 
et al. on a oxygen concentration of a;=0.35 was performed 
on a sample with a much higher Tc than discussed here, 
implying a larger hole doping in the Cu02 planesj^ It 
follows that the samples studied here have a hole concen- 
tration much closer to the critical value of Pc=0.055 than 
previous studies. 

The quasielastic scattering is similar to that in the 
monolayer La2-a;Sra;Cu04 near the critical doping of 
Pc=0.055 and so illustrates that the properties are not 
unique to the monolayer cuprates. One difference is 
that we observe commensurate magnetic peaks whereas 
the monolayer cuprates display incommensurate peaks. 
Since the resolution was adequate to detect a similar elas- 
tic incommensurate wave vector, it may be that the in- 
commensurability in underdoped YBCOe.ss is smaller or 
suffers intrinsic broadening because of short oxygen chain 
order. Nonetheless we have shown that the elastic corre- 
lation lengths and quasielastic moment are similar to the 
monolayer cuprate system. The total spectral weight at 
energies below 30 meV energy transfer has a similar mag- 
nitude to that in the insulating cuprates. We find that as 
the critical hole concentration is approached from above, 
slowly dynamic antiferromagnetic short-range order de- 
velops, as in a glass, at the expense of superconductivity. 

A spectrum with two energy scales similar to that 
measured here in hole-doped YBCOe.as also occurs in 
electron-doped Pro.88LaCeo.i2Cu04_5 (PLCCO) with 
Tc=21 K.^^'*^" This adds credence to our suggestion that 
the two energy scales with central mode as observed in 
the YBCOg+a;, LSCO, and PLCCO systems are a generic 
feature and common to underdoped high-temperature su- 
perconductors. 

The low-energy excitation spectrum in YBCOg.ss is 
very different from that of more heavily doped samples. 
The low-energy inelastic scattering is commensurate with 
no sign of a well-defined resonance peak. In the Ortho-II 
YBCO 6. 5 superconductor the incommensurate modula- 
tion is dynamic within the hour-glasSi-§ A suppression of 
magnetic scattering below ^ 16 meV below Tc, indicates 
the formation of a superconducting spin gap. No such 
gap or suppression is observed in YBCOg.ss, but instead 
an inexorable growth through Tc. 

Based on previous doping studies, we conclude that 
the lack of a clear resonance and spin-gap cannot be at- 
tributed to the presence of structural disorder. It may 
be possible that disorder present in the chains could lo- 
calize charge and have an effect on transport properties 
in the Cu02 planes. Effects of structural disorder on 
the superconductivity have been noted in excess oxygen 
La2Cu04+j, by Lee et al^ Doping impurities into the 
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YBCOe+a; superconductor has also been shown to in- 
duce long-ranged magnetic order. For cobalt doping, the 
cobalt ions enter the chains introducing not only struc- 
tural disorder but also antiferromagnetic order, even in 
heavily hole doped samples The order was long-ranged 
characterized by sharp Bragg peaks when scanned along 
both the [H,H,0] and [0,0, L] directions. Long-range order 
is patently not present in our YBCOe.as samples which 
display short correlation lengths both along c and within 
the a — b planes and a highly organized 3D spin pattern 
that indicates that significant structural disorder is ab- 
sent. While impurity and disorder effects have been char- 
acterized in heavily doped YBCO to have subtle effects 
on the lineshape and gap structure, they do not dramati- 
cally alter the lineshape and certainly do not remove the 
resonance and the spin-gap. Therefore, disorder cannot 
account for the large differences we observe in YBCOg.as 
and more heavily doped YBCOg.s. 

While the doped holes do not have a strong effect 
on integrated strength of low-energy excitations, they 
nearly destroy the high-energy excitations above ~ 200 
meV as they enter the pseudogap near the zone bound- 
aryi^ The low-energy fluctuations (characteristic of long- 
wavelength fluctuations) and the central peak are there- 
fore sensitive to spin disorder introduced through the 
doping of holes. It is only at high energies where the 
decay of of spin excitations into particle hole pairs is 
energetically possible, that the spin excitations become 
sensitive to electronic degrees of freedom. We note that 
the spin properties measured here are very different from 
those of heavily doped samples. This observation may be 
in accord with thermal conductivity experiments which 
have found a metal insulator transition in the YBCOe+a; 
system for x near 0.5i^ Indeed, a detailed study of the 
superconducting spin gap as a function of hole doping in 
YBCO suggests that the superconducting gap decreases 
rapidly for small hole doping 



from the excitations. We find overall that the behavior 
of a cuprate superconductor close to its critical doping 
is dramatically different from that of the heavily doped 
superconductors and involves a large shift of the spectral 
weight to low and quasielastic energies. 
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IX. APPENDIX: ABSOLUTE CALIBRATION 
USING BRAGG PEAKS 

To obtain an estimate of the strength of the quasielas- 
tic magnetic moment in YBCOg.as we carried out, in 
addition to the phonon calibration, an absolute calibra- 
tion against the nuclear Bragg peaks. This estimate al- 
lows a direct comparison with ordered antiferromagnets 
to determine how much of the Cu^'*' moment is ordered 
quasi-statically. Statically here means including the slow 
quasielastic peak since there is no strictly elastic scatter- 
ing in YBCOg.ss. The Bragg peak structure factor was 
determined from rocking scans in crystal angle using 
the relation 



VII. CONCLUSIONS 

We have presented a detailed study of the static and 
dynamic magnetic properties of YBCOe.as (Tc=18 K). 
Two distinct timescales exist, one a fast relaxational ex- 
citation with an energy linewidth of 2r=5 meV, and a 
second slowly dynamic centra mode that appears static 
at low temperatures on the neutron timescale. In con- 
trast to what is observed in more heavily doped samples, 
we do not observe a clear signature of superconductivity 
in the temperature dependence of the dynamic suscepti- 
bility, indicating a decoupling of the magnetic Cu^+ spins 
from the paired quasiparticles but not from the holes. For 
temperatures greater than 50 K, the magnetic scattering 
is well described by w/T scaling with the spectral energy 
width scaling linearly with temperature. A breakdown 
of scaling is observed at low temperatures and energies 
where the central peak develops and the spins gradually 
freeze into a spin glass state removing spectral weight 



/ # Imeas.W = ^7;^. (9) 

J sm[20) 

To carry out the calibration we plot the calculated struc- 
ture factor of the nuclear Bragg peaks against the mea- 
sured structure factors to obtain a calibration constant. 
The resulting curve is plotted in Fig. llSl and results from 
a fit to the following equation to correct for secondary 
extinction 

I 1 2 -^ca/cu/aied I 

\^ measured] Z ; fa' \^^) 

i -I- P\r calculated] 

The fitted calibration constant is a and (3 denotes a fac- 
tor to account for extinction. In the limit of small /3, 
Eqn. [To] agrees with the expression of Shamoto et al^ 
Given the complicated geometry of our sample mount 
consisting of seven co-aligned crystals, this simple model 
fits remarkably well. 



20 



lOr 



CO 



YBCO 



6.35' 



^ Bragg Peak Calibration Curve 
8 

/ 

5 
4 



3 



Ik 



2 I 
I 
I 

k 



.250 

calculated 



100 150 
(Arbitrary Units) 



200 



FIG. 18: The measured nuclear structure lactor squared as a 
function of calculated structure factor squared is plotted for 
YBC06.35. The data points are based on the integrated in- 
tensity calculated from transverse scans through Bragg peaks. 
The curve is the result of a fit to the empirical formula in the 
text that accounts for secondary extinction. 



traction method used throughout this study to obtain 
the elastic correlated magnetic scattering. By compar- 
ing the rocking scan through the (1/2, 1/2, 5) (shown in 
Fig. [TO)) position with the Bragg peak intensity, we ob- 
tain an ordered square moment of /i^=0.010±0.002 /i^ 
in YBC06.35 compared with ^ 0.47 /i^ measured by 
Shamoto et al. in insulating YBCOg.is. Therefore, the 
square of the static moment has decreased by an order 
of magnitude from insulating YBCOg.is to superconduct- 
ing YBCOe.ss- A similar decrease in the ordered moment 
has been observed in the single layer LSCO system as a 
function of hole dopingjii 
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The structure factor for the integrated intensity for two 
copper spins per unit cell is given by 

|FmP = ^^mV/^(Q)s^(Q). (11) 

where (7ro)^/4 is equal to 73 mbarns, g is the Lande 
factor set to 2, i?^(Q) is the bilayer structure factor pre- 
viously defined, and fj, is the value of the ordered mo- 
ment. The factor of 2/3 comes from the polarization fac- 
tor for isotropic spins as found from our polarized neutron 
analysis. The integral of the rocking scan of the central 
peak was obtained by subtracting from a 9 K scan a 
non-magnetic background at 80 K. This is the same sub- 



FIG. 19: A transverse {tp) scan through the correlated mag- 
netic scattering at (1/2, 1/2, 5). The data was obtained by 
subtracting a high temperature background taken at 80 K. 



The method used here to calibrate the elastic scatter- 
ing assume the peak is resolution limited along the lon- 
gitudinal direction. This is clearly not the case and may 
account for the slight underestimate of the elastic mo- 
ment when compared with the calibration analysis based 
on an acoustic phonon. Nevertheless, both methods agree 
well and provide justification for the quantitative analysis 
presented here. The value for the total ordered moment 
also provides a useful comparison with the insulator and 
other cuprate systems. 
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